EXPLOSIVE  DRIVEN  FERROELECTRIC  GENERATORS 


Ya.  Tkach 

Institute  of  Electromagnetic  Research 
Kharkov,  Ukraine 

sit  @  iemr.vl.kharkov.ua 


S.  Shkuratov  ,  J.  Dickens,  M.  Kristiansen 

Department  of  Electrical  Engineering 
Texas  Tech  University 
Lubbock,  TX 

sshkuratov@ppl.ee.ttu.edu 

L.L.  Altgilbers 

Advanced  Technology  Directorate 
U.S.  Army  Space  and  Missile  Defense  Command 
P.O.  Box  1500 
Huntsville,  AL  35807 

Larry.  Altgilbers  @  smdc.army.mil 


P.T.  Tracy 

Tracy  Physical  Sciences  Research 
Huntsville,  AL 

Ptracyl458@aol.com 


Abstract 

Explosive  Driven  Ferroelectric  Generators 
(EDFEGs)  are  compact  power  sources  that  have 
been  considered  for  use  as  seed  sources  for 
magnetocumulative  generators,  as  well  as  prime 
power  sources  [1,2].  Shock  waves  generated  by 
high  explosives  are  used  to  shock  depolarize 
ferroelectric  materials,  which  results  in  a  voltage 
pulse  being  delivered  to  a  load.  These  generators 
have  been  experimentally  investigated  at  Texas 
Tech  University.  Data  from  these  experiments 
was  used  to  benchmark  a  code  developed  at  the 
Institute  of  Electromagnetic  Research.  In  this 
paper,  a  description  of  the  simulation  and  a 
comparison  of  the  experimental  and  simulation 
results  will  be  presented. 

I.  DESCRIPTION  OF  MODEL 

As  the  shock  wave  passes  through  the  polarized 
module,  made  of  ferroelectric  ceramics  (in 
particular  lead-zirconate  titanate  (PZT)),  its 
volume  is  divided  into  two  zones  (Fig.  1) 
differing  by  such  parameters  as  bulk 
polarization,  permittivity,  and  conductance. 

These  zones  are  referred  to  as  the  “compressed” 
and  "uncompressed”  zones,  where  the 
compressed  zone  is  that  through  which  the  shock 


wave  has  already  passed  and  the  uncompressed 
zone  is  that  through  which  the  shock  wave  has 
not  passed.  The  equivalent  circuit  diagram  for 
the  longitudinal  EDFEG  is  presented  in  Fig.  1, 
where  Ci  and  C2  are  the  capacitances  of  the 
uncompressed  and  compressed  zones, 
respectively. 

In  building  the  model,  the  following 
assumptions  were  made: 

•  There  is  a  single  planar  shock  wave. 

•  Bulk  compression  of  the  ferroelectric 
(or  to  be  more  precise  ferroceramic) 
material  was  not  taken  into  account, 
since  experimental  results  indicate  that 
this  value  does  not  exceed  0.05. 

•  Shock  polarization  inherently  bears  an 
inertial  characteristic;  that  is,  domain 
rearrangement  is  a  kinetic  process. 

•  The  ferroelectric  material  is  a  linear 
dielectric  in  both  the  compressed  and 
uncompressed  zones. 

Taking  these  assumptions  into  consideration,  the 
general  set  of  ordinary  differential  equations  for 

the  current  in  the  reactive  load  /,  [A]  and  the 
charges  on  the  load  capacitor  Qn  [Cl]  and  PZT 
module  Ql  [Cl]  in  the  EDFEG  is 
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where  /?  (Ohm),  L  (H)  and  C/(F)  are  the 
resistance,  inductance,  and  capacitance  of  the 
load. 

The  capacitance  of  the  ferroelectric 
module  C(t)  is  described  by  using  the  standard 
model  for  a  layered  parallel  plate  capacitor, 
where  one  layer  is  the  compressed  zone  and  the 


other  layer  is  the  uncompressed  zone.  For  this 
type  of  capacitor,  the  distance  between  the 
plates,  l,  must  be  less  than  the  radius  of  the 

plates,  Sl,27U  .  Therefore,  the  capacitance  of 
the  ferroelectric  module  shown  in  Fig.  1  is: 


Cit) 


£l)£\ 
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where  S  (m2)  is  the  area  of  the  end  plates  of  the 
ferroceramic  module,  Et  is  the  permittivity  of  the 
uncompressed  zone,  e2 (p)  is  the  permittivity  of 
the  compressed  zone,  80  is  the  permittivity  of 
free  space,  /  is  the  total  length  of  the  ferroelectric 
module,  r  is  the  radius  of  the  module,  and  Vs  is 
the  velocity  of  the  shock  wave  in  the 
ferroelectric  material,  which,  in  the  absence  of 
substantial  changes  in  the  state  of  the  material, 
can  be  assumed  to  be  equal  to  the  velocity  of 
sound  in  the  material.  The  permittivity  of  the 
shock-compressed  ferroceramic  is  a  complex 
function  of  the  pressure  in  the  shock  wave. 

The  total  electric  charge  released  at  the 
end  plates  of  the  ferroceramic  module  during  the 
time  it  takes  for  the  shock  wave  to  travel  through 

the  module  is  Q  =  (j(p)S  .  Since  the  free 

charge  surface  density,  <J(p)  ,  released  at  the 
end  plates  is  equal  to  the  difference  in  the 
specific  bulk  polarization  in  the  compressed  and 
uncompressed  regions  of  the  material. 


for  /  <  Vst; 


O  =  Pl~  P2 ,  the  amount  of  charge  released  at 

any  given  moment  in  time,  Q(t(t)  ,  is 

proportional  to  the  depolarized  volume  of  the 
module  and  is  described  by  the  expression: 

00(t)  =  <T(p)vstr1. 

Taking  the  derivative  with  respect  to  time  yields 
the  depolarization  current: 

d<j(p)Vstll 


I0(t)  =  dW- 


clt 


where  0(t )  takes  into  account  the  “switching 
on”  of  the  polarization  current  at  the  origin, 
which  corresponds  to  the  moment  at  which  the 
shock  wave  enters  the  ferroceramic  module,  and 
the  “switching  off’  at  the  moment  the  shock 

wave  exits  the  module  ( tf  =  IV s  * )  and  is 
approximated  by  the  expression: 


(3) 

(4) 
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which  is  convenient  for  numerical  calculations, 
since  it  is  a  smooth  function.  Since  the 


calculation  of  t  ,(p)  using  kinetic  theory  is 


) 


1  +  tanh 


(5) 


the  end  plates  leaks  through  the  compressed 
region  of  the  module  forming  a  leak  current, 

I  leak  ( P)  ■  I11  the  case  of  a  longitudinally  driving 


relative  complex,  experimentally  measured 
values  of  0.05  -  0.4  ps  are  used  [3,4].  These 
values  decrease  significantly  as  the  pressure  in 
the  shock  wave  increases. 

The  electric  conductance  of 
ferroceramics  sharply  increases  under  shock 
compression  and  part  of  the  charge  released  at 


force,  the  leak  takes  place  over  the  entire  surface 
of  the  module  and.  consequently,  does  not 
depend  on  time.  The  leak  current  can  be  found 
from  the  field  strength  in  the  shock-compressed 
region  of  the  module  by  using  Ohm’s  law: 


Jleak(Q\ )  —  '^2(P)E2',  I  leak  ( Q\  )  —  Jleak(Q\)S  ~  ^2(P)Q\  _^2^0  ] 


(6) 


where  X2  (Ohm  m)  is  the  specific  conductance  of 
the  compressed  ferroceramic  material.  The 
conductance  increases,  as  the  pressure  in  the 
shock  wave  increases,  due  to  the  increase  in  the 
number  of  free  carriers  because  of  electron 
tunneling,  ionization,  and  other  phenomena. 

The  set  of  ODE  in  Eq.  (1 )  describes  the 
operation  of  a  longitudinally  driven  ferroceramic 
module,  where  the  shock  wave  velocity  and 
spontaneous  polarization  vectors  are  either 
parallel  or  anti-parallel,  before  the  onset  of  bulk 
breakdown  in  the  compressed  region  of  the 
module.  Generally  speaking,  breakdown  can 
occur  in  both  the  compressed  and  uncompressed 
regions,  but  it  starts  first  in  the  compressed 
region  because  the  electric  strength  of  the 
ferroceramics  is  less  than  that  in  the 
uncompressed  region.  This  is  related  to  the 
formation  of  local  breakdown  areas  in  the 
compressed  zone  due  to  the  impact  of  the  shock 
wave  on  powder  grain  boundaries,  defects. 


dislocations,  and  air-filled  cavities  generated 
during  the  baking  process.  Thus,  Eq.  ( 1)  is 
restricted  to  the  time  domain  prior  to  the  start  of 
bulk  breakdown  in  the  compressed  zone: 

Ebr>Q1(t)[Se2(p)e0]~1;  ii) 

where  Eh  [kV/m]  is  the  electric  strength  of  the 

ferroceramic  material  compressed  by  the  shock 
wave. 

Since  the  set  of  equations  in  Eq.  ( 1 )  is 
stiff  and  cannot  be  efficiently  solved  with  the 
required  precision,  they  are  normalized  by 
introducing  the  reduced  variables: 

T  =  t/t\  q  =  Q/Q',  and  i  =  It  / Q;  (8) 

where  t  =  IV s  '  and  Q  —  <y(p)  S .  After 

making  the  appropriate  substitutions,  Eq.  ( 1) 
becomes: 
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where: 
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^  leak  (<7i)  =  {£2(p)£oY'  Mp)hv 
The  normalized  bulk  breakdown  condition  is: 


Ebr><T(p)q1(T)[£2(p)£ o]"1. 


II.  RESULTS 


This  set  of  equations  was  solved  numerically 
using  the  Gear  and  Bulrich-Stoer  numerical 
methods  with  a  relative  error  for  both  methods  of 

less  than  10  5  for  all  variables. 


Substituting  in  the  parameters  for  the 
EDFEG  provided  by  Texas  Tech  University,  the 
model  was  used  to  calculate  the  output  voltage  of 
their  EDFEGs.  Since  some  of  the  parameters, 

such  as  T  . ,  (7,  £-, ,  and  X2  ,  are  not  known  at 
shock  pressures,  they  were  corrected  based  on 


earlier  shots  and  used  in  calculations  for  later 
shots.  It  should  be  noted  that  the  actual  values  of 
these  parameters,  after  corrections,  are 
reasonable  based  on  the  values  measured  for 
PZT  ceramic  materials  [3].  A  comparison  of  the 
calculated  to  the  experimental  data  for  two  shots 
is  presented  in  Fig.  2.  Experiments  show  that 
shock-compressed  ferroceramic  modules  can 
generate  pulses  with  amplitudes  up  to  8  kV  in 
the  resistive  part  of  the  load  and  energies  per  unit 
volume  of  module  of  0. 1  -  0.4  J/cm3.  The 
duration  of  the  pulse  depends  on  the  shock 
transit  time  through  the  module,  the  “skewness” 
of  the  shock  front,  and  the  relaxation  time  of  the 
ferroceramic  material.  Those  factors,  which 
probably  limit  the  output  energy  and  the 
amplitude  of  the  output  pulse,  are  the  electric 
breakdown  strength,  which  is  approximately  3 
kV/mm,  of  the  ferroceramic  material  and  the 
leakage  current  passing  through  the  compressed 
portion  of  the  module. 
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Fig.  1  Diagram  illustrating  the  propagation  of  a  shock  wave  through  a  longitudinally  driven  ferroceramic 
module  and  field  distributions  in  the  module. 


Figure  2.  Comparison  of  Experimental  to  Simulated  Output  Voltages  of  the  EDFEG  for  Two  Different 
generators  tested  at  TTU. 


